Centrifugal compression of soft particle packings — theory and experiment 
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An exact method is developed for computing the height of an elastic medium subjected to centrifu- 
gal compression, for arbitrary constitutive relation between stress and strain. Example solutions are 
obtained for power-law media and for cases where the stress diverges at a critical strain - for example 
as required by packings composed of deformable but incompressible particles. Experimental data 
are presented for the centrifugal compression of thermo-responsive N-isopropylacrylamide (NIPA) 
microgel beads in water. For small radial acceleration, the results are consistent with Hertzian 
elasticity, and are analyzed in terms of the Young elastic modulus of the bead material. For large 
radial acceleration, the sample compression asymptotes to a value corresponding to a space-filling 
particle volume fraction of unity. Therefore we conclude that the gel beads are incompressible, and 
deform without deswelling. In addition, we find that the Young elastic modulus of the particulate 
gel material scales with cross-link density raised to the power 3.3 ± 0.8, somewhat larger than the 
Flory expectation. 

PACS numbers: 83.80.Hj, 82.70.Gg, 62.20.D-, 46.65. +g 
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Interest in colloidal suspensions has been spurred by 
the advent of optical tools to image and manipulate 
behavior at the particle scale PJH]. While much re- 
search focusses on hard-sphere and charged systems, 
other work concerns the behavior of elastic microgel par- 
ticles composed of a swollen polymer network such as N- 
isopropylacrylamide (NIPA) [SHE]- Since such particles 
are soft and can swell/de-swell in response to variation 
of temperature, pH, or salt concentration, they are ideal 
as model systems for experiments at very high volume 
fractions where the particles are pressed together so that 
the material behaves as a jammed solid-like paste [9Hl6|. 

The mechanical behavior of a jammed packing ulti- 
mately originates in the elastic nature of the constituent 
particles. For microgel particles this depends on crosslink 
density and swelling state, and ought to form a crucial 
part of sample characterization. Similar considerations 
apply to packing of bubbles, cells, grains, etc. For large 
enough particles, the deformation of individual particles 
may be visualized under applied load - either of a packing 
[TThTU] or of a gel single bead [2"Ull23| . However, this is 
not feasible for submicron-scale colloidal microgel beads. 

In this paper, we demonstrate how the elasticity of a 
medium may be characterized by centrifugal compres- 
sion, and we illustrate our method with experiments on 
w 1 /im diameter NIPA microgel beads. Not surpris- 
ingly, theories for centrifugal compression have been pro- 
posed earlier. Ref. 24 approximates the overall sample 
compression in terms of an average pressure across the 
medium. Ref. |25j computes an approximate compres- 
sion profile for power-law media. By contrast we develop 
an exact prediction for sample height versus radial ac- 
celeration, for arbitrary stress versus strain constitutive 
relation. Our general theory applies to packings of small 
as well as large particles and, also, to any elastic medium 
such as an aggregated suspension [251 |57] that can be 



compressed either by centrifugation or gravity. Never- 
theless our primary interest here is in a non-cohesive 
random packing of spheres, for which we discuss how to 
relate bulk to particle-scale elasticity in light of recent 
advances regarding non-affine deformation. 



I. THEORY OF SAMPLE COMPRESSION 

Three main ingredients are required in order to pre- 
dict sample height H versus angular rotation speed ui for 
analysis of experimental data. This includes statements 
of mass conservation and force balance, as well as a con- 
stitutive model for compressive stress versus strain based 
on the elastic nature of the particle packing. Some of the 
key quantities for this task are defined on the schematic 
diagram of the experiment in Fig. [T] The sample itself 
consists of N particles in a suspending fluid of density pf . 
To encompass all possibilities of particle deformation and 
deswelling, we define m and v respectively as the mass 
and volume of particulate material in each particle. For 
example if the particles are incompressible homogeneous 
droplets, then m is the mass of the droplet and v is its 
volume. But if the particles are porous, then m and v 
do not include contributions from the suspending fluid 
within the pores. Therefore, in any case, m and v and 
the particulate material mass density p p — m/v are all 
constant and do not change no matter how the particles 
deform or deswell. 

In the limit of vanishing w, the particles are unstrained 
but close-packed at a volume fraction <j) c and a number 
density n c that are constant throughout the volume of 
the packing. The packing extends a radial distance or 
"height" H c inwards from the bottom of the container, 
with the supernatant fluid "above" - closer to the rota- 
tion axis. At nonzero u>, the particle packing compresses 




FIG. 1: (Color online) Schematic depiction of a soft particle 
packing under centrifugal compression. For rotation at angu- 
lar speed ui around the vertical axis as labeled, it comes to 
mechanical equilibrium at radial "height" H above the bot- 
tom of the sample. For this particular sample and centrifuge 
geometry, the cross sectional area A at depth z below the 
supernatant fluid-packing interface is constant except at the 
very bottom. 



to a smaller height H, and the strain 7 varies with the 
radial "depth" z below the supernatant fluid, so that ev- 
erywhere the elastic stress gradient counteracts the cen- 
trifugal buoyancy pressure gradient. In many centrifuges 
the sample tube is not perpendicular to the rotation axis, 
as in the depiction of Fig. [I] It is important to empha- 
size that the distance R between the rotation axis and 
the supernatant-packing interface, the "height" H, and 
the "depth" z, are all measured radially and not along 
the length of the tube. 



A. Mass conservation 

A thin slice of unstrained sample, with area A and 
thickness L, contains n c AL beads. If compressed to L — 
AL, the number density n increases but the number of 
beads nA(L — AL) is unchanged. Therefore the number 
density is 



1-7 



(1) 



where 7 = AL/ L is the compressive strain. If the entire 
sample has constant area A, and is compressed in ra- 
dial height from H c to H under centrifugation, then the 
total number of beads in the entire sample is similarly 
unchanged, n c AH c = J Q nAdz. The global expression 
of mass conservation is thus 



H r = 



H dz 

1-7' 



(2) 



More generally if the cross-sectional area varies across 
the sample, for example as at the bottom of the tube 
depicted in Fig. [Tl then the expression of mass conserva- 



sample tubes of very large length to diameter ratio, so 
that A can be considered as uniform and Eq. (pi) is accu- 
rate. Note that all these expressions of mass conservation 
hold whether sample compression is due to deformation 
of particles at constant volume, or deswelling of particles, 
or some combination of both. 



B. Force balance 

The particles in the thin slice of volume AAz depicted 
in Fig. [l] rotate with radius R + z and hence must ex- 
perience a net force that points radially inward and that 
equals total particle mass mnAAz times radial accelera- 
tion a r — lo 2 (R+z). Assuming that the sample boundary 
is frictionless, this net force arises both from the elastic- 
ity of the packing and from centrifugal buoyancy. The 
compressive strain 7 and elastic stress a zz increase with 
radial "depth" z, so the elastic force on the thin slice 
is A(da zz /dz)Az. The mass of fluid displaced by par- 
ticulate material within the slice is pfnvAAz, so the 
centrifugal buoyant force is given by Archimedes as this 
times the radial acceleration. Altogether, the statement 
of Newton's second law is 



,d<r 



(pfnvAAz)a r + A Az = (nmAAz)a 

dz 



(3) 



Here the right-hand side is mass of particles in the slice 
AAz times their acceleration, and the left-hand side is 
the sum of forces acting on the particles; the first force is 
due to the surrounding fluid and the second is due to the 
surrounding particles. Note that the factor AAz cancels, 
so that the only depth-dependent terms are the number 
density, given by Eq. (HI as n = n c /(l — 7), the elastic 
stress gradient da zz /dz, and the radius of motion (R + z) 
appearing in a r . Also, the constant n c (m — Pfv) is iden- 
tified as <j) c Ap, where <j) c is the volume fraction occupied 
by the particles at close-packing and Ap is the density 
difference between particle and fluid material. Thus the 
final expression for force balance simplifies to 



dcr* 



dz 



c Apu 2 (R + z) 

1^7 • 



(4) 



for any container shape and for any combination of par- 
ticle deformation and/or deswelling. 



C. Formal solution 

As a constitutive model we suppose that the elastic 
stress of the compressed particle packing may be written 
in general as 



Ya(rr) 



(5) 



r-H 



tion is V c = L [A/(l — 7)]dz. In our experiments, we use 



where Y is a materials property with dimensions of force 
per area and 5(7) is a dimensionless function of the com- 
pressive strain 7. We also suppose that sample height is 



small compared to the centrifuge radius, y < H ^ R. 
Then the force balance equation Q becomes 



(6) 

(7) 



Note that k is defined as a reciprocal length that char- 
acterizes both the materials and the radial acceleration. 
With a change of variables, the force balance equation 
may be integrated by parts: 



ds 
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dz 
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For a specific form of 5(7) characterizing the elastic na- 
ture of the packing, Eq. pj, the integration in Eq. ([9]) 
is to be performed and the result is to be expressed not 
in terms of 7 but ra ther in terms of s; this defines the 
function f(s) in Eq. (10) and its inverse f^ 1 (z). 



Finding the function f(s) defined by Eqs. (9 10) is tan- 
tamount to finding the height of the sample, if the cross- 
sectional area is constant. This can be seen by using the 
force balance equation (pi) to re-express the mass conser- 
vation equation (pi): 



H r = 



H ds 



1 



,1 1-7' 
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(11) 

(12) 
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Multiplying by k and taking the inverse gives the final 
formal result for radial packing height as a function of 
rotation speed: 



H = f(kH c ). 



(14) 



To recap, this solution assumes that the sample container 
has a constant cross-sectional area, and that the radius of 
circular motion is large compared to the packing height. 
However it makes no assumptions about whether sample 
compression is due to deformation or dcswelling of the 
particles. And it does not rely on explicit computation 
of strain versus radial depth. 



relations of potential experimental interest. The first and 
simplest is the elastic stress is a power-law of the strain, 



r 7 a , 



(15) 



so that the dimensionless function defined by Eq. (fsl) is 
s("/) = 7°. The value a = 1 corresponds to a linear 
spring-like medium and the value a = 3/2 corresponds to 
a packing of Hertzian spheres with small deformations. 
For the general case it is straightforward to carry out the 
integration in Eq. (p)]) and simplify to kz = 7° [1— 70/(1 + 
o)]. It is also straightforward to invert for 7 = s 1 / and 
re-express as kz — s[l — s 1 / a a/(l + a)], the right-hand 
side of which defines kf(s). According to Eq. (l4| ) the 
sample height is thus H = f(H c k) = (l/k)(Hck][l — 
(i/ c fc) 1 / a a/(l + a), which we express as 



H 



= 1 



(H c k) 



1/0 



(16) 



Since k ex oj 2 R 7 by definition in Eq. ([7]), the fractional 
decrease in sample height, and also the average strain 
of the entire sample (7) = 1 — H/H c , vary with the 
experimentally-accessible control parameters as a linear 
function of (uj 2 RH c ) 1 / cl and a proportionality constant 
that depends on materials parameters. Though Eq. ( 16 ) 



is remarkably simple, we emphasize that it is an exact 
solution, good for any radial acceleration or amount of 
sample compression. 

The form of Eqs. ( 14|16 ), and also of the predictions 
of sample height in the following sections for different 
constitutive laws, suggest that data may be conveniently 
analyzed in terms of a plot of H/H c versus the length 
x = uj 2 RH c J g. This should cause data collapse for differ- 
ent initial sample heights, and thus serves as a nice check. 
Then fits may be made to H/H c = 1- [a/(l + a)](Ka:) 1 / Q 
where the fitting parameter k is an inverse length de- 
fined so that H c k = kx. According to this definition and 
Eq. ([7| , the elastic constant of the bulk medium is then 



Y = (f> c Apg/n 



(17) 



The quality of the fit is bound to be good for small com- 
pression, and to give a reliable value for Y . Since Eq. ( fl6| ) 
is an exact solution, any deviation of the fit from the data 
at larger compression is directly related to a deviation of 
the actual constitutive relation from a zz = Y-f a . This 
makes centrifugal compression a sensitive probe of the 
elastic nature of the packing. For example, the sample 
height for an actual sample will surely not decrease to 
zero as predicted by Eq. ( 16 ) for strong but finite cen- 



trifugation, due to stiffening at large strains. 



D. Power-law elasticity 



E. Linear with maximum strain 



In this and the following two sub-sections, we use the 
above formalism to predict sample height versus radial 
acceleration for three specific stress-strain constitutive 



There must exist a maximum strain 7 m at which all 
fluid is expelled and the packing is pure particulate ma- 
terial. If the particulate material is incompressible, then 



the elastic stress must diverge at this maximum strain. 
A simple constitutive law that is linear at small strain 
and that diverges at j m is 



1 - l/l-n 



(18) 



Though ad-hoc, this form has the dual virtue of being in- 
vertible for strain versus stress and of being integrable in 
Eq. pi. Taking advantage of these features, the resulting 
packing height is computed exactly to be 



H 

H~, 



= 1 - 7, 



7m 



In 



1 



H c k 

7m 



(19) 



where the reciprocal length k is defined by Eq. Q, as 
before. For gentle centrifugation, H c k/^/ m <C 1, this may 
be expanded as H/H c = 1 - {H c k)/2 + (H c k) 2 /{3j m ) 
plus higher order terms; note that the leading behav- 
ior is identical to Eq. (16) for a = 1. For strong cen- 
trifugation, H c k/~f m 3> 1, the average strain of the 
entire sample approaches 7 m as H/H c = 1 — j m + 
~f m \n[H c k/j m }/(H c k/j m ) + 7 m /(i? c fc/7 m ) 2 plus higher 
order terms. This prediction may be a reasonable expec- 
tation for emulsions, for which the droplets are incom- 
pressible and are often assumed to interact as repulsive 
linear springs. 



F. Hertzian with maximum strain 

For a packing of solid spherical particles that are elastic 
but incompressible, we reason as above that the stress 
must be Hertzian at low strains and must diverge at some 
finite maximum strain r y m . The simplest such form we 
can conceive, that is also analytically tractable within 
the context of the formal solution for packing height vs 
compression, is 



Y 1 



3/2 



l-(7/7m) 3 / 2 ' 



(20) 



As above, this expression may be inverted for strain ver- 
sus stress, and it also may be integrated in Eq. Q. The 
resulting packing height is computed exactly to be 



H 



= i - 



7n 



7) 



5/2 



TT/VS 



W 3H c k 
2V3 tan" 1 (1 + 2W)/VS) 

-ln[(l + W + W 2 )/(l-W) 2 ] J, (21) 



W = 



H r k 



3/2 + Hr k 



1/3 



(22) 



The reciprocal length k is defined by Eq. ([7]), as be- 
fore. For gentle centrifugation, H c k/j m <C 1, the sample 
height prediction behaves as H/H c = l-(3/5)(H c k) 2 / 3 + 
(if c fc) 5 / 3 /(47 m 3 / 2 ) plus higher order terms; note that the 



leading behavior is identical to Eq. (JT6J) for a = 3/2. For 
strong centrifugation, it may also be verified that the 
sample compression asymptotes to H/H c = 1 — j c , as 
expected. 



G. Sphere and packing elasticities 

In this last subsection on theory, we discuss the connec- 
tion of microscopic particle properties to the macroscopic 
stress-strain relations for the case of elastic spheres. In 
particular, how does the measured value of the parameter 
Y defined by a zz — F7 3 / 2 depend upon the elasticity of 
the sphere material? To begin we recall the classic calcu- 
lation by Hertz for two elastic spheres of equal diameter 
d, brought into contact such that their centers are a dis- 
tance h closer together than the sum of their radii. For 
small deformation, the repulsive force is 



Ed 2 



3(1 



3/2 



(23) 



where E is the Young elastic modulus and v is the Pois- 
son ratio of the sphere material. A derivation of this 
expression is given by Landau and Lifshitz |28j . where it 
culminates in their Eq. (9.15). It is also given by Walton 
[29], where it culminates in his Eqs. (2.5 & 2.15), ex- 
pressed as combinations of Lame parameters that reduce 
to B = (1 - v 2 )I{-kE) and C = v{\ + u)/(irE). 

The main contribution in Walton's paper [5S] is com- 
putation of the effective elastic moduli of a random pack- 
ing of non-cohesive elastic spheres, assuming affine de- 
formation and averaging over a fixed set of randomly- 
oriented contacts at which the actual stress and strain 
fields are computed. The pressure required to achieve 
uniform compression 7^ = 7 is given in his Eq. (3.19) as 
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3/2 



3tt(1 



ZE 
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3/2 



(24) 
(25) 



where 4> is the volume fraction of spheres and Z is the 
average number of contacts per sphere. For "perfectly 
smooth [sic]" frictionless spheres that support no shear 
traction across the contact area, the stress for uniaxial 
compression j zz — 7 is given in Walton's Eq. (3.26) as 



2 .r" = Y-Y 



4>ZE .. ,, _ 
67r(l — v 
ZE 
67^^(1-^)' 



3/2 
) 3/2 . 



(26) 

(27) 



For "infinitely rough \sicY spheres that support any 
amount of shear traction across the contact area, the elas- 
tic stress is larger by a factor of (3 — 2f)/(2 — v). Walton 
also finds the anisotropy of the stress: cr xx / 'a zz = 1/4 for 
smooth spheres and *V(12 — 81/) for rough spheres. 



The connection between particle and packing behavior 
was broadly explored by O'Hern et al. |30j using numer- 
ical simulation. There the particles are frictionless disks 
or spheres that repel with central force F — (e/d)(h/d) a , 
where e is a microscopic energy scale, d is the particle 
diameter, and h is the compression. Different a, different 
dimensionality, and different particle size distributions 
are all examined. The pressure P, shear modulus G, 
bulk modulus K, and coordination number Z are always 
found to scale with volume fraction as 



p = 


-- P o ((f>-0 c ) a 


(28) 


G = 


-- Go^-^T- 1 ' 2 


(29) 


K = 


= Ko^-^y- 1 


(30) 


Z = 


-- z c + z (4> - cp.) 1 ' 2 


(31) 



Here K = acj)P follows from K = -VdP/dV = 
4>dP/d4>, Z c equals twice the dimensionality of the sys- 
tem, and 4> c — 0.639 ± 0.001 for monodispcrse three di- 
mensional spheres. Some of these scaling relations were 
found in previous [31H33] and subsequent [3H [35] simula- 
tions for special cases. An important general conclusion 
from all these studies is that the microscopic deformation 
is nonaffine - the particle positions adjust away from the 
macroscopic deformation field in order to reduce individ- 
ual compression. This causes the scaling exponent of the 
shear modulus to be +1/2 larger than the affine expec- 
tation, so that resistance to shear is dramatically lower 
near <j) c . By contrast the exponent for the bulk modulus 
is unchanged, but K is typically 10-30% smaller than 
if the motion were affine. These findings suggest that 
Walton's prediction, Eq. (26), is an overestimate. 



It is instructive to directly compare the pressure pre- 
diction of Walton with the simulations of O'Hern, ct 
al. For monodisperse Hertzian elastic spheres in 3- 
dimensions, a = 3/2, <p c = 0.64, and Z c — 6, the simula- 
tion results quoted in Table I of [30] are P — 0.35(e/d 3 ), 
G = 0.1A(e/d 3 ), and Z — 7.7. Agreement of the sim- 
ulation force law with the Hertz Eq. (j23h requires that 



the energy scale be taken as e = Ed 3 /[ 



% Alto- 



gether this gives the expectation based on the O'Hern et 
al. simulation as P = Q.lYJ[E/{\ - v 2 )]{(f> - C ) 3/2 - As 
expected due to nonaffine motion, this is slightly smaller 
(by a factor of 0.76) than Walton prediction Eq. ( [25] ) 
evaluated with Z = 6 and <j> c = 0.64. Therefore, if the 
particles adjust under uniaxial compression so that the 
stress becomes isotropic, a lower bound on a zz would be 
about 0.76 times Walton's pressure prediction; this corre- 
sponds to 0.29 times his Eq. (26) prediction for Y. In the 



absence of further guidance, we will simply use Eq. ( 26 ) 



to deduce E from measurements of Y; the actual Young 
modulus could be up to three times larger. 



II. EXPERIMENT 

In this section we use the above theory to analyze 
the elasticity of a suspension of N-isopropylacrylamide 
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FIG. 2: (Color online) Particle volume V versus tempera- 
ture T, determined by dynamic light scattering from a dilute 
sample. The line is an empirical fit, V = (2.93 /im 3 )[l — 
T/(39.6°C)]. Above 35° C the volume collapses to 0.092 /jm. 
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FIG. 3: (Color online) Radial packing height H versus radial 
acceleration for samples at T = 20° C but with different initial 
heights, H c , as indicated by the points along the y-axis. The 
curves represent fits to Eq. ( 21 1 with H c held fixed. The 



statistical uncertainty in measuring H is about 0.1 mm, much 
less than the symbol size. 



(NIPA) microgel beads. The particles are synthesized 
via free-radical polymerization pj H31 [TSl I3"6"H3"5] . Briefly, 
NIPA monomer (Acros) and methylene-bis(acrylamide) 
crosslinker (Polyscience, Inc.) are mixed in aqueous so- 
lution. Ammonium pcrsulfate (Fisher) is then injected 
into solution to initiate the polymerization. Spherical 
particles are formed within an hour. Next these are thor- 
oughly washed and redistributed in a 1 mM sodium do- 
decylsulfate solution. In Ref. [16] we reported on mi- 
crofluidic measurements of the shear rheology for dense 
suspensions of the very same samples. 

For our samples the particle diameter is of order 1 
micron in diameter. But more importantly, as shown 
by the dynamic light scattering (DLS) results plotted in 
Fig. [2] the volume can decrease by more than a factor 
of 2 with a modest temperature increase. Over the tern- 
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FIG. 4: (Color online) Packing height versus acceleration, 
with both axes scaled so as to cause collapse, for samples at 
different temperatures, as labeled. Different color symbols 
represent different initial sample heights as in Fig. [3] The 
pink dotted lines are the initial decays, which are linear on 
such a plot for Hertzian particles. The solid curves are fits 
to Eq. (pi]), and asymptote to 0.64 ± 0.01 as indicated by the 
horizontal dashed line. 
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FIG. 5: (Color online) Temperature dependence of the param- 
eters k and maximum strain -y m , obtained from fits of Eq. ( 21 1 
to normalized compression data as demonstrated in Fig. 141 
The assumed stress-strain relation, Eq. (20 1, is Hertzian at 
low strains, a zz — Y7 3 ' 2 , and the sample resists compression 
beyond a maximum strain y m . The sample elasticity scales 
as Y = 4> c £i.pg/ k. The 7 m results are constant to within 
uncertainty, and average to 0.36 ± 0.01 as indicated by the 
horizontal dashed line. The error bars are set by the degree 
of collapse and accuracy of the fits in Fig. HI 



perature range studied here, the particle volume exhibits 
an approximately linear dependence with temperature as 
noted in the figure caption. At high temperatures, above 
approximately 35° C, the particles suddenly collapse to a 
temperature-independent volume of 0.092 /im. For hard 
colloidal spheres sterically stabilized by graphted poly- 
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FIG. 6: (Color online) Young's elastic modulus E of the par- 
ticulate material plotted versus (a) temperature and (b) par- 
ticle volume. The values of E are deduced from the fitting 
parameters k in Fig. [5k using Y — 4> c Apg/K and Walton's 
relation in Eq. ( |26[ ). The dashed line in (b) is a power-law, 
with slope as labeled. It translates to the dashed curve in (a) 
using the empirical fit to particle volume versus temperature 
in Fig. [| 



mer, the "hydrodynamic radius" given by DLS is some- 
what larger than the hard-sphere radius given by electron 
microscopy. Here, the NIPA particles are sterically sta- 
bilized by dangling chains that emerge from crosslinking 
sites within the gel. It is not possible to remove the 
beads from solution and measure their size by electron 
microscopy. However, since the rheological features re- 
ported in Ref. [16! exhibit power-law behavior in (f> — (f> c 
where 4> c = 0.635 ± 0.003 coincides well with the volume 
fraction of randomly close-packed spheres, we conclude 
that the hydrodynamic radius given by DLS corresponds 
closely with the actual physical radius of the particles. 

As the particles shrink with increasing temperature, 
the density of crosslinking sites necessarily increases and 
the elastic modulus hence increases in some way that 
we wish to determine. Our experimental protocol is as 
follows. We load different volumes of a stock suspen- 
sion of particles into six 50 fiL glass capillary tubes (ID 
= 0.8 mm, length = 10 cm), and seal at both ends 
with optical glue. The number density of particles is 
0.284//im 3 , as determined by counting particles in a 
three-dimensional confocal microscopy image of a sample 
with known dilution. Loading is done at room tempera- 
ture, where the particles occupy a volume fraction of ap- 
proximately 40%. The sample volumes are chosen so that 
once the particles sediment to close packing, the range of 
initial radial heights H c (see Fig. II) spans between about 
3 and 6 cm. This is long enough that the shape of the 
bottom of the tubes plays no role. We load all the tubes 
into a thermostated centrifuge (Marathon 21000R), and 
let the particles settle at a fixed angular rotation speed uj. 
This centrifuge has a rotation radius of R = 10 cm and 
angle away from vertical of 30°; the temperature range 
is about 10-25°C and is held constant to ±0.1°C. To 
determine the time needed for complete settling to me- 
chanical equilibrium, we measure the height versus time 
for a range of rotation rates by periodically removing the 
tubes and tracking the sediment-supernatant interfaces. 
The settling equilibrates within a few hours for high ro- 



tation rates, but takes up to several days for low rotation 
rates. Once the samples have settled to their equilibrium 
heights we place the tubes in a holder on an optical ta- 
ble, and photograph using a digital Nikon D70 camera. 
We analyze the images for the length of the sediment 
and convert to the radial height, H, via the appropriate 
trigonometric factor (See Fig. IT]). The statistical uncer- 
tainty is about AH = 0.1 mm. Then we agitate the 
samples to redistribute the particles, and repeat at a dif- 
ferent rotation rate. It should be noted that the packing 
eventually rebounds elastically to random close packing, 
but this process occurs on a scale of many hours to days 
while our height measurements take only minutes. Wc 
also note that the particles clearly return to their origi- 
nal spherical shape with no plastic deformation: we can 
repeat the experiment with a redistributed sample or a 
fresh sample and produce the same result. 

Example data for radial height versus angular acceler- 
ation are shown in Fig. [3] for all six tubes at temperature 
T = 20° C. The points at uj = are not from centrifuge 
measurements, but instead are the expectations for H c 
based on particle number density, particle size, sample 
volumes, and a random close packing volume fraction of 
4> c = 0.64. Since the colloidal gel particles are presum- 
ably Hertzian, and since the height data do not decrease 
toward zero, we fit to Eq. (21 1 keeping H c fixed. These 



fits are all good, and interpolate smoothly between the 
expected H c at w = and the centrifuge results at u > 0. 
This gives confidence in our sample characterization and 
centrifugation measurement procedures. 

We now attempt to collapse the compression data ac- 
cording to the general expectation in the theory section 
for any constitutive law. In particular, in Fig. [4] we 
plot the normalized radial height H/H c versus the length 
x = uj 2 RH c /g raised to the 2/3 power, for all six tubes 
and for three different temperatures. As such, the y-axis 
represents the observation while the x-axis consists of a 
combination of the three control parameters that could 
be varied in experiment. We see in Fig. HI that this does 
indeed cause good collapse of the height data at each 
temperature. Also as expected, at higher temperatures 
the particles are smaller and hence stiffer and less com- 
pressed. Furthermore the initial decay is linear on such 
a plot, in accord with Hertzian behavior at small strains. 

Next we fit the collapsed data, for all tubes at a given 
temperature, to the prediction Eq. (21 ) for a sphere pack- 



ing that is Hertzian at small strains and that cannot 
be strained beyond some maximum j m . The first fit- 
ting parameter is the reciprocal length, n, defined as dis- 
cussed earlier so that H c k = nx in Eq. (21); the value 



of k is determined by the elasticity of the medium and 
sets the slope of the initial linear decay seen in Fig. [4] 
The second fitting parameter is the maximum strain 7 m , 
which sets the asymptotic value of the scaled height as 
H a /H c = 1 — 7„j at high rotation speed. The fits at all 
temperature are as satisfactory as those shown in Fig. [4] 
Therefore the elasticity of the packing is adequately de- 



parameters have the intended physical meaning. 

The fitting parameters k and j m are plotted versus 
temperature in Fig. [5j The top plot shows how k de- 
creases with temperature, and is analyzed in the follow- 
ing paragraph. The bottom plot shows that the maxi- 
mum strain is independent of temperature to within ex- 
perimental uncertainty, and averages to r y m = 0.36±0.01. 
This gives an asymptotic relative height of H a /H c — 
1 — 7m = 0.64 as shown by the dashed line in Fig. |4J Tak- 
ing the random close packing fraction as <j) c — 0.64 and 
assuming the particles deform at fixed volume without 
deswelling, the volume fraction (j> m at maximum strain 
would be given by <j) c H c = (j) m H a as <j) m — 1. This is 
space filling, and implies than any compression H / H c 
below the dashed line at 1 — -f m = 0.64 in Fig. [4] could 
only be accomplished by deswelling. Since the actual 
data at all temperatures are always above this limit, and 
appear to approach it smoothly from above, we conclude 
that the particles do not deswell and are effectively in- 
compressible. This is consistent with reports that the 
Poisson ratio of bulk samples of swollen NIPA is close to 
1/2 [3SJ |30]. It is also consistent with the observation 
that ps 10 MPa of applied pressure is needed for notice- 
able deswelling of NIPA microgels [31], while the pres- 
sure here does not exceed 0.03 MPa. Stated differently, 
we estimate that a radial acceleration of about 10 6 m/s 2 
(200 times our maximum) would be required to induce 
deswelling and, thus, to cut off the divergence assumed 
in Eq. ([20|. See Ref. [22l [42] for the deswelling of poly- 
electrolyte gels under compression. 

We now analyze the k results in Fig. [5]l for the elas- 
ticity of the NIPA particles. Recall that the Hertzian 
elastic constant in a zz = F7 3 / 2 is given by Eq. (17) 
as Y = (pcApg/n. Here the density difference Ap be- 
tween particles and water at T — 20° C is found to 
be 0.08 g/cm 3 by measuring the terminal sedimentation 
speed of single spheres at room temperature, and equat- 
ing gravity to Stokes forces. Density values for different 
temperatures are then deduced using the particle volume 
versus temperature data of Fig. [2] After converting n to 
Y, we finally deduce the Young elastic modulus accord- 
ing to Eq. {26} as E = 6tt(1 - v 2 )Y/{(f)Z), using v = 1/2 



scribed by the empirical form of Eq. ( 20 1 , and the fitting 



and evaluating the denominator at random-close pack- 
ing, 4> — 4> c = 0.64 and Z = 6. The results are plotted 
in Fig. p] versus (a) temperature and (b) particle volume. 
As expected, E increases with temperature since the par- 
ticles shrink. The order of magnitude is tens of kPa, as 
found previously for bulk NIPA gel samples 39, 40, 43 

US]. 

The scaling of E vs V plotted in Fig.|6]3 is fit well by a 
power-law of E ~ i/y3.3±0.8 Since the number of cross 
links in the gel beads does not change as they shrink with 
temperature, the same exponent x — 3.3 ± 0.8 holds for 
the scaling of E with cross-link density. For this Flory 
theory applied to neutral gel [46H48] predicts a value of 
x = 2.25 in good solvent and x — 3 in poor (theta) sol- 
vent. Experimental values mentioned in |48j are x = 2.3 
and 2.4 for good solvents and x = 3.0 and 3.7 for theta 



solvents. In our case, as the temperature is increased, not 
only cross-link density increases but also, solvent qual- 
ity decreases [39] . which leads to an increased exponent 
value. Furthermore, Flory theory supposes that the gel 
network deswells in an affine way. Affine deformation can 
be impaired two-fold: heterogeneities in the cross-links 
distribution may arise from the synthesis mechanisms 
[3Hllin]i while entanglements between cross-links may be 
created as the network shrinks. Departure from affine 
behavior results in an increase in the E versus cross-link 
density exponent beyond the x — 2.25 value [J5]. Exper- 
imental evidence for such a non-affine shrinking of NIPA 
gels were reported in |40) where a harder skin is found 
to form at their surface upon temperature increase. In 
our case, such a skin would have to be thicker than the 
scale of deformation under compressions in the Hertzian 
regime. And finally, our samples are not perfectly neu- 
tral as assumed for the Flory value of X = 2.25; rather, 
the NIPA beads have a very slight negative charge, and 
the surfactant is anionic. This situation can be compared 
with polyelectrolyte gels swollen in brine, where Hertzian 
behavior and x = 3.5 was reported in Ref. ;21 from com- 
pression of individual beads. Altogether, the decrease of 
solvent quality and the possibilities of non-affine shrink- 
age and of charging effects are all consistent with our 
experimental measure of x — 3.3 ± 0.8 for the scaling 
exponent of E with cross-link density. 



III. CONCLUSION 



is greatest at the bottom of the sample and decreases 
to zero at the top. We hope that the example solutions 
developed for power-law media, and for harmonic and 
Hertzian media with a maximum strain, will be of use to 
experimentalists. These same predictions hold for media 
that are so soft that gravity causes measurable compres- 
sion, by replacing ui 2 R with g. 

To illustrate, and to characterize the elastic nature of 
particles in a suspension of thermoresponsive gel beads 
of interest for shear rheology experiments |16j . we pre- 
sented a series of measurements of sample height versus 
rotation speed, filling height, and temperature. The re- 
sults are in good agreement with Eq. (21 1, the exact so- 



lution for a Hertzian medium with maximum strain as 



specified by the constitutive law of Eq. (20 1. The data 



thus are analyzed in terms of the linear elastic Young 
modulus E of the gel material, and in terms of the max- 
imum strain j m beyond which the medium cannot be 
compressed. The former is demonstrated to vary as a 
large power of cross-link density, and was used in Ref. |16j 
to non-dimensionalize shear rheology data. The latter is 
a constant, j m — 0.36 ± 0.01, whose value implies that 
the gel beads are incompressible and do not deswell when 
deformed; this conclusion is important in Ref. [16] for 
knowing the packing fraction above jamming. Valuable 
particle-scale information such as this is now straightfor- 
ward to obtain by centrifugal compression measurements, 
analyzed using the theory presented here. 



Here we presented both theory and experiment for the 
compression of an elastic medium under centrifugation. 
The formalism culminates in Eq. ( 14 ) for the total sam- 
ple height expressed in terms of a function f(s) defined 
by Eqs. (9p0). Remarkably, this represents an exact so- 



lution for arbitrary stress-strain elastic constitutive law 
and was found without having to first solve explicitly for 
the height dependence of the compressive strain, which 



Acknowledgments 

We thank P. E. Arratia, F. Lequeux, A.J. Liu, and A. 
G. Yodh for helpful conversations; we thank A. Alsayed, 
A. Basu, and Z. Zhang for helping synthesize the parti- 
cles. This work was supported by the National Science 
Foundation through grants MRSEC/DMR05-20020 and 
DMR-0704147. 



A. van Blaaderen and P. Wiltzius, Science 270, 1177 
(1995). 

J. C. Crocker and D. G. Grier, Journal of Colloid and 

Interface Science 179, 298 (1996). 

D. G. Grier, Nature 424, 810 (2003). 

V. Prasad, D. Semwogerere, and E. R. Weeks, Journal of 

Physics-Condensed Matter 19, (2007). 

H. Senff and W. Richtering, Journal of Chemical Physics 

111, 1705 (1999). 

B. R. Saunders and B. Vincent, Advances in Colloid and 
Interface Science 80, 1 (1999). 

R. Pelton, Advances in Colloid and Interface Science 85, 

1 (2000). 

D. M. Heyes and A. C. Branka, Soft Matter 5, 2681 

(2009). 

M. Cloitre, R. Borrega, and L. Leibler, Physical Review 



Letters 85, 4819 (2000). 
[10] M. Cloitre, R. Borrega, F. Monti, and L. Leibler, Physical 

Review Letters 90, 068303 (2003). 
[11] J. R. Seth, M. Cloitre, and R. T. Bonnecaze, Journal of 

Rheology 50, 353 (2006). 
[12] D. A. Sessoms, I. Bischofberger, L. Cipelletti, and 

V. Trappe, Philosophical Transactions of the Royal Soci- 
ety A - Mathematical Physical and Engineering Sciences 

367, 5013 (2009). 
[13] J. Mattsson, H. M. Wyss, A. Fernandez-Nieves, 

K. Miyazaki, Z. B. Hu, D. R. Reichman, and D. A. Weitz, 

Nature 462, 83 (2009). 
[14] Z. X. Zhang, N. Xu, D. T. N. Chen, P. Yunker, A. M. 

Alsayed, K. B. Aptowicz, P. Habdas, A. J. Liu, S. R. 

Nagel, and A. G Yodh, Nature 459, 230 (2009). 
[15] K. Chen, W. G Ellenbroek, Z. X. Zhang, D. T. N. Chen, 



P. J. Yunker, S. Henkes, C. Brito, O. Dauchot, W. van [33 

Saarloos, A. J. Liu, et al., Physical Review Letters 105, 

025501 (2010). [34 

[16] K. N. Nordstrom, E. V. Vcrneuil, P. E. Arratia, A. Basu, 

Z. Zhang, A. G. Yodh, J. P. Gollub, and D. J. Durian, [35 

arXiv:1007.4466vl (2010). 
[17] T. Lachhab and C. Weill, European Physical Journal B [36 

9, 59 (1999). 
[18] T. S. Majmudar and R. P. Behringer, Nature 435, 1079 [37 

(2005). 
[19] M. M. Bandi, M. K. Rivera, F. Krzakala, and R. E. Ecke, [38 

arXiv:0910.3008vl (2009). 
[20] D. C. Andrei, B. J. Briscoe, P. F. Luckham, and D. R. [39 

Williams, Journal de Chimie Physique et de Physico- [40 

Chimie Biologique 93, 960 (1996). 
[21] A. Knaebel, S. R. Rebre, and F. Lequeux, Polymer Gels [41 

and Networks 5, 107 (1997). 
[22] A. Knaebel and F. Lequeux, Polymer Gels and Networks 

5, 577 (1997). [42 

[23] R. D. Egholm, S. F. Christensen, and P. Szabo, Journal 

of Applied Polymer Science 102, 3037 (2006). [43 

[24] R. Buscall, Colloids and Surfaces 5, 269 (1982). 
[25] W. H. Shih, W. Y. Shih, S. I. Kim, and I. A. Aksay, [44 

Journal of the American Ceramic Society 77, 540 (1994). 
[26] K. T. Miller, R. M. Melant, and C. F. Zukoski, Journal 

of the American Ceramic Society 79, 2545 (1996). [45 

[27] J. J. Lietor-Santos, C. Kim, P. Lu, A. Fernandez-Nieves, 

and D. Weitz, European Physical Journal E 28, 159 [46 

(2009). 
[28] L. D. Landau and E. M. Lifshitz, Theory of Elasticity, [47 

vol. 7 of Course of Theoretical Physics (Butterworth- 

Heinemann, Oxford, 1986), 3rd ed. [48 

[29] K. Walton, J. Mech. Phys. Solids 35, 213 (1987). 
[30] C. S. O'Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, [49 

Phys. Rev. E 68, 011306 (2003). 
[31] D. J. Durian, Phys. Rev. Lett. 75, 4780 (1995). 
[32] D. J. Durian, Phys. Rev. E 55, 1739 (1997). 



H. A. Makse, N. Gland, D. L. Johnson, and L. M. 
Schwartz, Phys. Rev. Lett. 83, 5070 (1999). 
H. A. Makse, N. Gland, D. L. Johnson, and L. Schwartz, 
Phys. Rev. E 70, 061302 (2004). 

W. G. Ellenbroek, E. Somfai, M. van Hecke, and W. van 
Saarloos, Phys. Rev. Lett. 97, 258001 (2006). 
A. M. Alsayed, M. F. Islam, J. Zhang, P. J. Collings, and 
A. G. Yodh, Science 309, 1207 (2005). 
M. Reufer, P. Diaz-Leyva, I. Lynch, and F. Scheffold, 
European Physical Journal E 28, 165 (2009). 
P. Yunker, Z. X. Zhang, and A. G. Yodh, Phys. Rev. 
Lett. 104, 015701 (2010). 
J. Hirotsu, J. Chem. Phys. 94, 3949 (1991). 
W. C. Lin, K. R. Shull, C. Y. Hui, and Y. Y. Lin, J. 
Chem. Phys. 127, 094906 (2007). 

J. J. Lietor-Santos, B. Sierra-Martin, R. Vavrin, Z. B. Hu, 
U. Gasser, and A. Fernandez-Nieves, Macromolecules 42, 
6225 (2009). 

S. Vervoort, S. Patlazhan, J. Weyts, and T. Budtova, 
Polymer 46, 121 (2005). 

T. Takigawa, T. Yamawaki, K. Takahashi, and T. Ma- 
suda, Polymer Gels and Networks 5, 585 (1997). 
T. R. Matzelle, D. A. Ivanov, D. Landwehr, L. A. Hein- 
rich, C. Herkt-Bruns, R. Reichelt, and N. Kruse, J. Phys. 
Chem. B 106, 2861 (2002). 

T. R. Matzelle, G. Geuskens, and N. Kruse, Macro- 
molecules 36, 2926 (2003). 

P. G. de Gennes, Scaling Concepts in Polymer Physics 
(Cornell University Press, NY, 1979). 
M. Rubinstein and R. Colby, Polymer Physics (Oxford 
University Press, NY, 2003). 

S. P. Obukhov, M. Rubinstein, and R. H. Colby, Macro- 
molecules 27, 3191 (1994). 

G. Nisato, F. Schosseler, and S. J. Candau, Polymer Gels 
and Networks 4, 481 (1996). 



